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normalization constant, equation (A12) 

Airy function, equation (All) 
speed of light, fm/s 

matrix element of qq pair interaction, equation (1) 
terms in approximation to Airy function, equation (A28) 
Planck’s constant, MeV-s 
total meson spin 

strength of potential, equations (2) and (7) 

orbital meson momentum 

associated Laguerre polynomial 

orbital quantum number 

mass of the eigenstate, equation (1) 

magnetic quantum number 

mass of quark or antiquark ( q may be b, c, or t) 

transformation variables, equations (B22) 

radial quantum number 

radial node number, equation (4) 

associated Legendre polynomial 

radial component of wave function, fm"3/2 
zeros of Airy function 
radial distance, fm 
intrinsic meson spin 

reduced radial component of wave function, fm -1 /^ 
constant additive potential, equation (1) 
potential function, equations (2) and (7) 
spherical harmonics 
arbitrary function 
gamma function 

zenith component of wave function 
transformation variables, equations (A6) and (A8) 
reduced mass 

azimuthal component of wave function 
complete wave function 
angular frequency 



Introduction 


With the present interest in understanding quark 
confinement, a large amount of theoretical work 
based on nonrelativistic models (ref. 1), relativistic 
models (ref. 2), and field theoretic models (ref. 3) 
has been done to obtain particle mass spectra of 
mesons as bound states of quarks. According to 
the quark model, mesons are composed of a tightly 
bound quark-antiquark pair ( qq ) which interacts 
with an external probe like a single collective ob- 
ject. Of particular interest are the nonrelativistic 
quark models, which have proven to be very suc- 
cessful for describing the spectroscopy, decay am- 
plitudes, static properties, and binding forces of 
mesons and baryons (refs. 4 through 11). These non- 
relativistic models are most successful when they are 
applied to systems where the energy level spacings 
are much smaller than the energy equivalent of the 
constituent masses (i.e., for systems which are, in 
fact, nonrelativistic) such as mesons consisting of the 
heavy charm, bottom, and top quarks. Through 
the use of nonrelativistic models, several indepen- 
dent lines of arguments suggest that the confinement 
of quarks may be based on an interaction which in- 
creases for large quark separation in proportion to 
the separation itself. This argument is assumed to be 
the dominating mechanism which explains the energy 
eigenvalues in the qq pair system for mesons. 

In this paper, the picture that charm, bottom, 
and top q and q are interacting with each other 
through a linear potential is used to analyze the 
particle spectrum of observed (charm, bottom) and 
predicted (top) mesons. Based on this picture, the 
quarks in the qq pair initially freely recede from each 
other in accordance with the ideas of asymptotic free- 
dom; however, as the separation distance increases, a 
linear potential emerges and a bound state is formed. 
It is then assumed that, as long as the quarks are in 
an energy eigenstate, their lifetime is reasonably long 
and their eigenstates can be completely characterized 
by a linear potential. 

Through the use of experimentally observed res- 
onances at 3096 and 3686 MeV/c 2 for the charmed 
meson, the two parameters of the linear potential k 
and Vo (the strength of the interaction potential and 
an additive constant) are determined. Then, in the 
center-of-mass system of the bound qq pair, the non- 
relativistic Schrodinger equation is exactly ( t — 0) 
and asymptotically (£ ^ 0) solved for a linear po- 
tential to obtain the remaining energy eigenvalues 
for charm and all eigenvalues of the bottom and top 
mesons. One then adds twice the quark mass plus a 
constant to these energy levels to obtain the meson 


mass as 

M n c 2 = 2m q c 2 + Vq + E n (1) 

where M n c 2 is the meson mass, m q c 2 is the con- 
stituent quark mass, Vo is a constant additive 
term, and E n is the matrix element of the qq pair 
interaction. 

The quark confinement potential has been treated 
for many years as a harmonic potential mainly be- 
cause of its closed-form solution in both the non- 
relativistic and relativistic regimes (refs. 12 through 
16). Here the exact nonrelativistic solutions for the 
eigenvalues and eigenfunctions of the harmonic po- 
tential are presented and the resulting energy levels 
are compared with those of the linear potential to 
validate the use of these two potential forms. 

Potential Models 

Linear Potential 

The three-dimensional linear potential is defined 
as 


v(r ) = kr (2) 

where k is the constant describing the strength of the 
interaction potential and r is the qq pair separation. 
With h in units of MeV-s and c in units of fm/s, 
the energy eigenvalues for an S state {l = 0) are ob- 
tained from the exact solution of the nonrelativistic 
Schrodinger equation as 


Efir — 


k 2 (hc)* 

pc 2 


-|l/3 


Rn 


where 


( 3 ) 


n r = n — 1 (n = l, 2, 3, ...) (4) 


and Rn represents the zeros of the Airy function 
along the negative abscissa. Appendix A describes 
the detailed derivation of equation (3) and its corre- 
sponding wave function. For i ^ 0, an exact solution 
does not exist; however, asymptotic solutions have 
been obtained to predict the energy levels. Among 
these, for arbitrary n and l. is the Cornwall equation 
(ref. 17) 
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and the Sheehy and Von Baeyer equation (ref. 18) Results 
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Note that in equation (5), E n i depends only on 
one combination of n and £. Therefore, there is a de- 
generacy between the states (n, t) and (n — 1, t + 2). 
In equation (6), E nr i with n r defined in equation (4) 
is derived from a modification to the Bohr quantiza- 
tion of circular orbit with the added radial excitation 
being treated as a harmonic approximation and by 
the use of half-integer angular momentum quantum 
number which, when applied to the confined heavy 
mesons, yields a reasonably accurate spectrum. 


Harmonic Potential 

The three-dimensional harmonic potential is de- 
fined as 


v{r) = ~kr 2 (7) 

where k and r are defined in equation (2). The energy 
eigenvalues for arbitrary n and £ are given by 

E nr t ~ (% n r + i + he ^^2 ^ (®) 

Appendix B describes the detailed derivation of equa- 
tion (8) and its corresponding wave function. 

Spin Consideration 

Since quarks have intrinsic spin of 1/2, the fol- 
lowing states are possible: , 3 5 q , l P\ ,Po ?P\ ?P2 > 

etc., where the notation is 
(25+1 w ith S as the intrinsic spin of the me- 
son, L as the orbital momentum, and J as the total 
meson spin. The notation 5, P, and D correspond 
to the L values of 0, 1, and 2, respectively. Note that 
all these states are observed; however, in this work 
the concentration is on S and P states as shown in 
figures 1 and 2, with the corresponding masses of 
states obtained from reference 19. Note that spin 
effects are* neglected, and average values for the 3 P 
states shown in the last column (figs. 1 and 2) are 
included. This column is actually not the average 
but rather the average of the mass of the 3 P2 state 
plus error and 3 Po state minus error, with the error 
quoted for the average 3 P state being just made to 
extend up to the mass of the 3 P 2 state plus error and 
down to the 3 Po state minus error. Thus the mass 
and the error assigned to the average 3 P state simply 
represent the range of the possible values of 3 Pq, 3 Pi, 
and 3 p 2 states. 


The usual procedure in quarkonium research is to 
treat the charm or bottom mass m c <? and m^c 2 as 
free parameters, and based on the particular poten- 
tial model in the nonrelativistic Schrodinger equa- 
tion, fit E n in equation (1) to two arbitrary experi- 
mentally observed levels to determine k and Vq. This 
means that there is a variety of possible fitting pro- 
cedures that one could try. Here the charm quark 
mass, at a typically accepted value of 1500 MeV/c 2 
(ref. 4), and the first two 3 Si state levels of char- 
monium at 3096 and 3686 MeV /c 2 are used in equa- 
tion (1) to obtain values of k and Vjp With the first 
two levels fitted to the experimentally observed lev- 
els, k and Vq are then used to predict other levels. 
With the linear potential model of equation (2), the 
calculated values for k and Vo are & = 1211 MeV/fm 
and Vq = 690 MeV. 

The results of this iterative procedure are shown 
in figures 3 and 4 and are compared with experi- 
mentally observed levels. Also included in figures 3 
and 4 are the results of the harmonic potential model 
(eq. (7)), which were calculated by the authors in an 
earlier paper (ref. 16). In these figures, levels ob- 
tained by the linear potential model are labeled with 
an L, whereas levels obtained from the harmonic po- 
tential model are labeled with an H. Note that, for 
the linear potential model, good agreement is found 
in both the charmonium and bottomonium spectra 
for S and P levels. In almost all cases, the linear po- 
tential model produces a much more accurate spec- 
trum than the harmonic potential model. The main 
advantage of the harmonic potential, however, is that 
the solution of the nonrelativistic Schrodinger equa- 
tion is much more easily obtained than for either the 
linear or more exotic potential models. Therefore, 
the harmonic potential model has traditionally been 
used to obtain initial estimates for energy levels as a 
first try. 

Recently, claims have been made that the sixth 
quark flavor, called top, has been found (refs. 20 
and 21) with a mass in the range of 30000 to 
50 000 MeV/c 2 . Choosing the top mass to be 
40000 MeV/c 2 , calculations for toponium levels 
based on the linear potential model are presented 
in figure 5 along with predictions from the harmonic 
potential model. 

As stated earlier, for i ^ 0, the nonrelativistic 
Schrodinger equation, with a linear potential, can- 
not be solved exactly. Thus asymptotic solutions of 
equations (5) and (6) were presented as good approx- 
imations for calculation of energy levels for arbitrary 
£ values. Table 1 presents the experimental results, 
exact results of equation (3), and approximate results 
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of equations (5) and (6) for charm, bottom, and top 
mesons, for the reader to draw conclusions on the 
applicability of equations (5) and (6). As is apparent 
from the table, in almost all cases, the Cornwall equa- 
tion (eq. (5)) produces a more accurate spectrum 
than the Sheehy and Von Baeyer equation (eq. (6)). 
This is partially because equation (6) was derived 
based on the modification to the classic Bohr quan- 
tization of circular orbit, with the added radial ex- 
citation being treated in a harmonic approximation 
(perturbation), and equation (5) was derived based 
on purely asymptotic arguments without any added 
perturbative term and produces accurate energy lev- 
els even for toponium which is a highly debated topic 
in the present physics community. 

Figures 6 through 11 present the linear and 
the harmonic potential wave function U nr ^[r\ for 
charm, bottom, and top, in units of fm -1 / 2 for 
n r = 0, 1, 2, and 3, and t = 0 and 1, where 
the radial node n r was defined by equation (4). 
For the linear case with i = 0, the wave func- 
tion was derived from the exact solution of the non- 
relativistic Schrodinger equation (eq. (A32)). For 
t 7 ^ 0, a fourth-order Runge-Kutta scheme was used 
to numerically solve the nonrelativistic Schrodinger 
equation and obtain the energy levels ( 3 P) and wave 
function. For the harmonic case with arbitrary l, 
the wave function was derived from the exact so- 
lution of the nonrelativistic Schrodinger equation 
(eq. (B30)). In all figures, U Ur t[r\ is plotted for the 
range 0.01 < r < 4 fm. As is apparent, the wave func- 
tion is mass dependent and becomes larger in ampli- 
tude and narrower in width with increasing mass. 


Finally, note that even though this and other 
works (refs. 10 and 22) yield agreement between the- 
ory and experiment with a linear potential model, 
there are other potential forms (such as the logarith- 
mic form of Richardson (ref. 23)), which also give 
reasonable descriptions. This is because the main re- 
quirement is that the potential be confined. Good 
fits are also obtained by considering the form of the 
potential at intermediate distances (refs. 24 and 25), 
which represents a correction to the linear confine- 
ment form. 

Concluding Remarks 

Linear and harmonic potential models in conjunc- 
tion with the nonrelativistic Schrodinger equation 
were used to obtain particle mass spectra of mesons 
as bound states of quarks. The main emphasis was 
on the linear potential where an exact solution of the 
S-state (£ = 0) eigenvalues and eigenfunctions and 
an asymptotic solution for the higher order partial 
waves {l > 0) were obtained. Cornwall’s asymptotic 
solution very accurately predicted all the energy lev- 
els including recently claimed toponium. Exact so- 
lutions of eigenvalue and eigenfunction for harmonic 
potential with arbitrary partial wave states were also 
obtained and were compared with linear potential re- 
sults. For quark confinement, the linear potential is a 
much superior model, even though the solution of the 
nonrelativistic Schrodinger equation is more difficult 
to obtain. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
May 31, 1988 
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Appendix A 

Derivation of Linear Potential Energy Levels and Wave Function 

In this appendix, a detailed derivation of equation (3) and its corresponding wave 
function is presented. 

The three-dimensional linear potential is given by 

v(r ) = kr (Al) 

where k is the strength of potential and equation (Al) is subject to the limiting conditions 

= { £ (r < 0) } <A2) 

Because of the spherical harmonic nature of this potential, the angular dependence of the 
wave function can be given by the spherical harmonics; therefore, 


1 /,(r,0,<f>) = R[r]Y £ m (0,<t>) (A3) 

where ip is the complete wave function, and only the radial component R[r] must be 
calculated, which satisfies the nonrelativistic Schrodinger equation according to 


d 2 £{£ + 1) yE ukr 

o 1“ 


dr 2 


n 2 n 2 


fl[r] = 0 


(A4) 


An analytical solution of equation (A4) for £ ^ 0 does not exist. However, for t = 0, 
equation (A4) becomes 


d 2 \iE fikr 

dr 2 h 2 h 2 

By defining new variables £ and rj 

c __ VE 

* h 2 

ak 


equation (A5) becomes 


R[r ] = 0 (A5) 

(A6a) 
(A6b) 

(A7) 

r)r- £ = (A8) 

where /? is a constant. Using equation (A8), after some algebra, equation (A7) becomes 

(A9) 

Equation (A9) with the condition that /? 3 = rj 2 becomes 


-~g2 -Z + r ) r J R \ r ] = 0 


The term rjr — £ can be written as 


' d 2 d 3 \ 

d>?~^ X ) R[r] = 0 
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' d 2 \ 

W] = ° 


(A10) 


The solution to equation (A10) is the Airy function which is closely related to the Bessel 
function of fraction order and, in integral form, is expressed as 


4 


i f°° 
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/2[A] = A Ai{ A] 


(All) 
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where A, the normalization constant, is calculated later. With the boundary condition at 
r = 0, 


we find 


or 


and at boundary r = 0, 


we find 


or 


R{(3\\ = R{r}t - f]| r=0 = 0 

(A13) 

1 

0$ 

II 

(A14) 

mi = r 

(A15) 

1 

II 

1 

[( 

(A16) 

II 
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with rj = gives 



(A18) 

>-(*)"* 

(A 19) 


Therefore, eigenvalues of E are proportional to the roots of Airy function, and thus at 
r = 0, 


R 


-(*) 


1/3 
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= Ai 


-u>) 


1/3 


= 0 


or 


En r — 


h 2 k 2 \ l/z 

I \Rn\ (n r = n — 1 = 0, 1, 2, ...) 


A 4 


(A20) 


(A21) 


Equation (A21) gives the energy levels in terms of n r and R n (zeros of Airy function). We 
can rewrite this result in terms of our units as 
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En r = 


k^(hc)^ 


,1/3 
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Now, since 


PX = rjr - ^ - 


n 2 
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or 
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h* 




1/3 /Vfc\ 1/3 

A = -R n + I ^ * r 


which means Rn r [r] can be expressed as 

Rn r [r] = A At 

where the normalization constant A is 
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Therefore the wave function R Ur [r] is 


ii 1/3 / nk 

RnM = A * _i? " + (** 


0 
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(A24) 


(A25) 


(A26) 


(A27) 


In order to evaluate equation (A27) numerically, use reference 26 and the identity 


Ai[z] = cif[z] - c 2 g[z\ 


(A28) 


where 


/w= S 3t (i) 
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Defining the reduced wave function {/[r] as 


U[r) = 


R[r] 


(A30) 


Equation (A29) becomes 


U„ r \r] = 


{ 


Cl/ 


-R n + (pfc/h 2 ) ^ r\ — C2gr 


Rn + ( 2 ) *■] | ( 71 - 3 / 2 ^ 1 / 2 ^ ) 


1/3 


(A31) 


With h in units of MeV-s and c in units of fm/s, equation (A31) becomes 


{ci/ [--Rn + (/ic 2 fc/(M 2 ) 1/3 r| -C 2 ? [-i?n t (^c 2 fc/(ftc) 2 ) 1/3 rj| 
t/n r [r] 


2 T 1/3 

(A32) 


Equation (A32) in units of fm 1 / 2 along with equation (A22) describes the wave function 
and energy levels of linear potential with t = 0 completely. 
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Appendix B 

Derviation of Harmonic Potential Energy Levels and Wave Function 

In this* appendix, a detailed derivation of equation (8) and its corresponding wave function 
is presented. 

The three-dimensional harmonic potential is given by 

v{r) = i kr 2 (Bl) 

where 

k — fiu ; 2 (B2) 

with 


m Q mq m Q 

u = — — = 

m q + rriq 2 

being the qq pair reduced mass and 


~er 

being the angular frequency of the oscillation with k as the spring constant. In appendix A, 
it was shown that the complete wave function can be written as the product of a radial 
component R[r] and a spherical harmonics <j>) as 


iP(r,O,ct>) = Rlr)Y £ m (0,<f>) 


(B3) 


and only the radial component R [r] must be calculated which satisfies the nonrelativistic 
Schrodinger equation according to 


dr 2 


By defining new variables £ and r] 
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equation (B4) becomes 
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dr] 2 rj 2 h(j 

A solution for R[r]] in equation (B6) must be found subject to the following limiting conditions: 


A reasonable choice for R[r]\ can then be 


{i) o) 
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therefore, 


= r/ (+1 exp 


(-H 


Pin] 


(B9) 


Replacing R[rj\ in equation (B6) with the one in equation (B9), equation (B6) becomes 


d 2 -PM 

dr ] 2 


+ 2 




dP[n] 

dr] 


+ 




P[n] = o 


(Bio) 


By defining a new variable / as / = rj 2 and after some algebra, equation (BIO) becomes 

P{f] = 0 (Bll) 




E-W*» 


or 


dp f[\ 2 J J \ df + 2f[V + 2j hu. 


P[/]=0 


If 


_ u 

x = £+ - 
2 


y=2 


Lv 2) hu\ 


Then equation (B12) becomes 


* p[ n + ^ f) m-y Plf]= o 


d/ 2 /' 


df f 


which is the standard confluent hypergeometric equation, and 


(B12) 

(B13a) 

(B13b) 

(B14) 

(B15) 


Equation (B15) indicates that in equation (B9) when rj — ► oo and R[rj] — ► exp^— 

equation (B15) can be expanded in a finite number of terms. This is true when y in 
equation (B13b) is zero or a negative integer; therefore, 



(n r = 0, 1, 2, 3, ...) 


(B16) 


or 

E nr l — ^2n r + £ + hco (B17a) 

and with equation (B2), E Ur i becomes 

E nr t = ( 2 n r + i + he (4^ V2 (B17b) 

which is simply the energy levels in terms of n r and i and expressed in units of MeV. Now, 
equation (B9) can be written as 
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RnrlWi = Ar ) e+1 ex P 




-n r , £+ T} 2 


(B18) 


where A , the normalization constant, is calculated later. 

In equation (B15), P ^—n r , i + §; rj 2 J is related to the associated Laguerre polynomials by 

rt+1/ 2 m _ (r[nr+^+i]) 2 „ r . . 3 


-nr (/I = V 1 , r ,V' P 
n r ir y. + j 


- n r , £+ 2 ; / 


(B19) 


(B20) 


Therefore, equation (B18) becomes 

Rn r e = A v e+1 exp (~^ 2 ) L l n T 1,2 {r} 2 ] 
where L^ l ^ 2 [ri 2 ] is defined as 

,<+i/2, 2 , = v ( nr+<+ i) ! tfn 

m=o (n r - m)\ (m + £ + ! m! 

From a computational point of view, the sum in equation (B21) is cumbersome to calculate. It 
however can be recast by defining two new variables N and M as 


(B21) 


IV — Tlf + £ 


(B22a) 


M = m + £ (B22b) 

Therefore, n r = N — £, m = M — £, n r — m = {N~£) — ( M — £) = N — M, and equation (B21) 
in terms of £, M, and N can be written as 


M=N 


(N+jY 


2\M—t 


jr«+l/2r 2i _ V- kHOZ 

ihi (IV — M)\ (M + 5)! 

Now, by using the identity 


1\,_ (2fc + l)l 

2 2k+1 k\ 


H>= 


M W 


Equation (B23) becomes 


,2 \M-t 


.<+1/2, ,, _ "f." [(2/V + l)V(2^>/V!)] (-,2 ).- 

nr WJ [(IV — A/)!] [(2M + 1)!/(2 2A * +1 M!)] 


which simplifies to 


,£+1/2, 2 , = ^ [(2IV + 1)!] (2 2A ^ +1 M!) (~r? 2 ) M - g 

^ ^ [(IV-M)!][(2M + l)!](2 2 ^+ 1 iV!) (M — I!)! 

Therefore, by using the results of equations (B22) and (B26), equation (B20) becomes 


(B23) 


(B24) 


(B25) 


(B26) 


n , _ jUL+1 exp (-IA M f (2IV + l)!(2 2Ai+1 M!) 

^ AT] exp ^ 2 n ) [(7V-M)!][(2M + 1)!](2 2JV + 1 IV!) [M - 1)\ 1 ’ 


(2 N + 1)!(2 2M+1 M!) 


„2 \M-t 
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where A can be calculated according to 



Rldr] dr = 1 


(B28) 


and by using equations (B5), (B22), (B27), and ( B28 ), the radial wave function R nr t{r\ becomes 


Rn r t = 


2n! 


( ^ V 1/4 

\y/fic 2 kj \r[n + ^+|] 


1/2 


l y/vPk) 


e + 1 


exp 


-l 


_ l f he A r 2 

2 \y/Jic^k) 


M=N 


[(2N + 1)\}(2 2M+1 M\) (hc/y/ji^iej r 


*jC [(/V — M)!][(2M + 1)!](2 2JV+ 1/V!) 




(M - f)! 


(B29) 


In equation (B29), R Ut l is in units of fm 3 / 2 , and, finally by using equation (A30) with h in 
units of MeV-s and c in units of fm/s, equation (B29) becomes 



m [(2AT + 1)!](2 2M+1 M!) (hc/s/nc 2 id) r 2 

X Xl ~ M ) ! 1K 2M + 1)!](2 27V+1 JV!) (M - t)\ 


(B30) 


Equation (B30) in units of fm */ 2 along with equation (B17b) describes the wave function and 
energy levels of harmonic potential completely. 
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Table 1. Experimental and Theoretical Resonances for Linear Potential 


n 

n r 

l 

Experimental 

resonances, 

MeV/c 2 

Calculated resonances, 

MeV )c i 

Linear 
potential 
(eq. (3)) 

Sheehy and 
Von Baeyer 
(eq. (6)) 

Cornwall 
(eq. (5)) 

Charmonium 

i 

0 

0 

3 096 

3 095 

3 172 

3 091 

2 

1 

0 

3 686 

3 684 

4 098 

3 684 

3 

2 

0 

4 030 

4 166 

5 023 

4 167 

4 

3 

0 

4 415 

4 592 

5 948 

4 594 

1 

0 

1 

3 485 

3 441 

3 463 

3 408 

2 

1 

1 


3953 

4 105 

3 935 

3 

2 

1 


4 399 

4 747 

4 386 

4 

3 

1 


4 802 

5 388 

4 792 

Bottomonium 

1 

0 

0 
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9 633 

9 687 

9 632 

2 

1 

0 

10 023 

10 030 

10 311 

10 032 

3 

2 

0 

10 355 

10 355 

10 934 

10 358 
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0 

10 573 

10 643 

11 558 

10 645 
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0 
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9 868 

9 883 

9 846 

2 

1 

1 
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10 213 

10 316 

10 201 

3 

2 
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10 748 

10 505 

4 

3 

1 
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11 181 
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Toponium 
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Figure 1. Experimental energy level spectra of charmonium. Nmnbers on each level are exact values in MeV/c 2 . 
Values are taken from reference 19. 
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Figure 2. Experimental energy level spectra of bottomonium. Numbers on each level are exact values in 
MeV/c 2 . Values are taken from reference 19. 
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Figure 3. Calculations of linear and harmonic potential model levels versus experimental levels for charmonium 
spectrum. 
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Figure 4. Calculations of linear and harmonic potential model levels versus experimental levels for bottomonium 
spectrum. 
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Figure 5. Calculated toponium spectrum. 
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